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Abstract Neutrino-matter cross sections and interaction rates are central to the 
core-collapse supernova phenomenon and, very likely, to the viability of 
the explosion mechanism itself. In this paper, we describe the major 
neutrino scattering, absorption, and production processes that together 
influence the outcome of core collapse and the cooling of protoneutron 
stars. One focus is on energy redistribution and many-body physics, 
but our major goal is to provide a useful resource for those interested 
in supernova neutrino microphysics. 
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1. Introduction 

Supernova explosions are one major means by which elements are 
injected into the interstellar medium and, hence, into subsequent gen- 
erations of stars. Therefore, supernovae are central to the chemical 
evolution and progressive enrichment of the universe. Most supernova 
explosions are the outcome of the dynamical collapse of the core of a 
massive star as it dies. Collapse creates high temperatures (> 1 MeV) 
and densities (10 7 g cm -3 < p < 10 15 g cm -3 ) and produces ("after the 
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dust settles") either a neutron star or a black hole. Under such extreme 
thermodynamic conditions, neutrinos are produced in abundance. The 
mechanism of core-collapse supernovae is thought to depend upon the 
transfer of energy from the inner core to the outer mantle of the iron core 
of the massive star. Neutrinos seem to be the mediators of this energy 
transfer. Therefore, to fully understand core-collapse supernova explo- 
sions one must have a firm handle on the physics of neutrino production, 
absorption, and scattering. 

In this paper, we summarize the neutrino-matter cross sections and 
the neutrino production rates in the core-collapse context. Some of this 
discussion can already be found in Burrows (2001). We do not attempt 
to explain the hydrodynamics of supernova explosions, but do try to 
present the relevant neutrino processes that play a role. For the former, 
the reader is referred to Burrows (2000), Burrows, Hayes, & Fryxell 
(1995), Liebendorfer et al. (2001a), Liebendorfer et al. (2001b), Rampp 
& Janka (2002), Buras et al. (2003b), and Thompson, Burrows, & Pinto 
(2003). 

In §2, we present a physical derivation of stimulated absorption and 
then in §3 we summarize the basic neutrino-matter cross sections. In §4, 
we discuss the neutrino-electron scattering kernel, along with a simple 
treatment of the collision integral. In §5, we provide the relativisitic for- 
malism for inelastic neutrino-electron and neutrino-nucleon scattering 
processes and energy redistribution for non- interacting nucleons. This 
is followed in §6 with a discussion of the alternate, more powerful, for- 
malism for determining differential interaction rates and redistribution in 
the many-body context, namely that of dynamical structure factors. The 
role of strong and electromagnetic interactions between nucleons and lep- 
tons is explored, as well as collective excitations of the medium. Source 
terms for electron-positron annihilation (§7), neutrino-anti-neutrino an- 
nihilation (§8), and nucleon-nucleon bremsstrahlung (§9) cap off our 
review of the major processes of relevance in core-collapse simulations. 

2. Stimulated Absorption 

The concept of stimulated emission for photons is well understood 
and studied, but the corresponding concept of stimulated absorption for 
neutrinos is not so well appreciated. This may be because its simple 
origin in Fermi blocking and the Pauli exclusion principle in the context 
of net emission is not often explained. The net emission of a neutrino is 
simply the difference between the emissivity and the absorption of the 
medium: 



Jnet — f]v K a l v , 



(1) 
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where K a is the absorptive opacity and I u is the specific intensity. All ab- 
sorption processes involving fermions will be inhibited by Pauli blocking 
due to final-state occupancy. Hence, rj u in eq. (1) includes a blocking 
term, (1 — J r I/ )(Bruenn 1985). T v is the invariant distribution function 
for the neutrino, whether or not it is in chemical equilibrium. 

We can derive stimulated absorption using Fermi's Golden rule. For 
example, the net collision term for the process, v e n <-> e~p, is: 

r f d3 Pve f d3 Pn f d 3 p p f d 3 p e (s^imA 

^ e " P J (2v? J {2*fJ (2vr)3 ) 

x E(u e n <-> e~p) (2vr) 4 S 4 (p l/e + p n - p p - p e ) , (2) 

where p is a four-vector and 

Z(v e n <-> e~p) = T Ve T n {\ - T e )(\ - T v ) - T e T v {\ - T n ){\ - F Ve ) . (3) 

The final-state blocking terms in eq. (3) are manifest, in particular that 
for the v e neutrino. Algebraic manipulations convert E(z/ e n <-> e~p) in 
eq. (3) into: 



E(v e n^ep) = T n (\ - T e )(\ - T t 



q(l -Fv^) 3~v, 

-l — j u 

F n {l-F e ){l-F p ) 



where 



jpcq = [ e (£„ e -( Me -/i))/3 + ^-1 (5) 

is an equilibrium distribution function for the v e neutrino and it has 
been assumed that only the electron, proton, and neutron are in thermal 
equilibrium, ft is the difference between the neutron and the proton 
chemical potentials. Note that in T^t there is no explicit reference to a 
neutrino chemical potential, though of course in beta equilibrium it is 
equal to p, e — fi. There is no need to construct or refer to a neutrino 
chemical potential in neutrino transfer. 
We see that eq. (4) naturally leads to: 

Jnet = {B v ~ h) = <{B V - I v ) . (6) 

Of course, B u is the black body function for neutrinos. This expres- 
sion emphasizes the fact that C Uen ^ e - p and J ne t are the same entity. If 
neutrinos were bosons, we would have found a (1 + Jv q ) in the denomi- 
nator, but the form of eq. (6) in which l v is manifestly driven to B v , the 
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equilibrium intensity, would have been retained. From eqs. (4) and (6), 
we see that the stimulated absorption correction to n a is 1/(1 — J7> q ). 
By writing the collision term in the form of eq. (6), with K a corrected 
for stimulated absorption, we have a net source term that clearly drives 
I u to equilibrium. The timcscale is 1/ck*. Though the derivation of 
the stimulated absorption correction we have provided here is for the 
v e n <-> e~p process, this correction is quite general and applies to all 
neutrino absorption opacities. 

Kirchhoff's Law, expressing detailed balance, is: 

n a = r] u /B u or n* a = rf v /B v , (7) 

where n' v is not corrected for final-state neutrino blocking. Furthermore, 
the net emissivity can be written as the sum of its spontaneous and 
induced components: 

where + or — is used for bosons or fermions, respectively. Eq.(7) can be 
used to convert the absorption cross sections described in §3 into source 
terms. 



(8) 



3. Neutrino Cross Sections 

Neutrino-matter cross sections, both for scattering and for absorp- 
tion, play the central role in neutrino transport. The major processes 
are the super-allowed charged-current absorptions of v e and v e neu- 
trinos on free nucleons, neutral-current scattering off of free nucleons, 
alpha particles, and nuclei (Freedman 1974), neutrino-electron/positron 
scattering, neutrino-nucleus absorption, neutrino-neutrino scattering, 
neutrino-antineutrino absorption, and the inverses of various neutrino 
production processes such as nucleon-nucleon bremsstrahlung and the 
modified URCA process (y e + n + n — > e~ + p + n). Compared with 
photon-matter interactions, neutrino-matter interactions are relatively 
simple functions of incident neutrino energy. Resonances play little role 
and continuum processes dominate. Nice early summaries of the vari- 
ous neutrino cross sections of relevance in supernova theory are given in 
Tubbs & Schramm (1975) and in Bruenn (1985). In particular, Bruenn 
(1985) discusses in detail neutrino-electron scattering and neutrino- 
antineutrino processes using the full energy redistribution formalism. 
He also provides a practical approximation to the neutrino-nucleus ab- 
sorption cross section (Fuller, Fowler, & Newman 1982; Aufderheide et 
al. 1994), though better physics soon promises to improve the treatment 
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of these cross sections substantially (Langanke & Martinez-Pinedo 2003; 
Pruet & Fuller 2003). For a neutrino energy of ~10 MeV the ratio of the 
charged-current cross section to the i/ e -electron scattering cross section 
is ~100. However, neutrino-electron scattering does play a role, along 
with neutrino-nucleon scattering and nucleon-nucleon bremsstrahlung, 
in the energy equilibration of emergent neutrinos (Thompson, Bur- 
rows, & Horvath 2000). 

Below, we list and discuss many of the absorption and elastic scatter- 
ing cross sections one needs in detailed supernova calculations. In §4, 
§5, and §6, we provide some straightforward formulae that can be used 
to properly handle inelastic scattering. The set of these processes com- 
prises the essential microphysical package for the simulation of neutrino 
atmospheres and core-collapse supernovae. 

3.1 v e + n — > e~ + p: 

The cross section per baryon for v e neutrino absorption on free neu- 
trons is larger than that for any other process. Given the large abun- 
dance of free neutrons in protoneutron star atmospheres, this process is 
central to v e neutrino transport. A convenient reference neutrino cross 
section is a , given by 



4G%(m e c 2 ) 2 AA 9 

a n = — ^ V h - L705 x 10 cm > (9) 
7r(nc) 4 



where Gf is the Fermi weak coupling constant (~ 1.436 x 10 49 ergs 
cm -3 ). The total v e — n absorption cross section is then given by 



<„^i^)(^)V(^)Y'w 

(10) 

where g& is the axial-vector coupling constant (~ —1.23), A np = m n c 2 — 
nipC 2 = 1.29332 MeV, and for a collision in which the electron gets all 
of the kinetic energy e e - = e Ve + A np . Wm is the correction for weak 
magnetism and recoil (Vogel 1984; §5) and is approximately equal to 
(l + l.le u Jm n c 2 ). At e Ve = 20 MeV, this correction is only ~ 2.5%. We 
include it here for symmetry's sake, since the corresponding correction 
(Wm) f° r neutrino absorption on protons is (1 — 7.1sp e /m n c 2 ), which 
at 20 MeV is a large —15%. To calculate k*, <t" n must be multiplied 
by the stimulated absorption correction, 1/(1 — ^ ), and final-state 
blocking by the electrons and the protons a la eq. (4) must be included. 
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3.2 v e + p — > e + + n: 

The total v e — p absorption cross section is given by 



1 + ^9 a \ ( e ^e — ^np 



1 



2 



A n p 



1/2 



(11) 



where e e + = — A np and is the weak magnetism/recoil correction 
given in §3.1. Note that is as large as many other corrections and 
should not be ignored. To calculate k*, <Tp eP must also be corrected for 
stimulated absorption and final-state blocking. However, the sign of p e — 
\jl in the stimulated absorption correction for v e neutrinos is flipped, as is 
the sign of p, e in the positron blocking term. Hence, as a consequence of 
the severe electron lepton asymmetry in core-collapse supernovae, both 
coefficients are very close to one. Note that the v e + p — > e + + n process 
dominates the supernova neutrino signal in proton-rich underground 
neutrino telescopes on Earth, such as Super Kamiokande, SNO, and 
LVD, a fact that emphasizes the interesting complementarities between 
emission at the supernova and detection in Cerenkov and scintillation 
facilities. 

3.3 is e A <-► A'er 

From Bruenn (1985) the total v e — A absorption cross section, is ap- 
proximated by 



1 - 



m e c 2 

Ev + Q' 



1/2 

Wblock 



(12) 

where W bloc k = (1 - f e ~) e^-^-Q )^ Q> = M A , - M A + A ~ fi n - 
fip + A, A is the energy of the neutron I/5/2 state above the ground 
state and is taken to be 3 MeV (Fuller 1982), and the quantities N p (Z) 
and N n (N) are approximated by: N p (Z) = 0, Z - 20, and 8 for Z < 20, 
20 < Z < 28, and Z > 28, respectively, and N n (N) = 6, 40 - N, and 
for N < 34, 43 < N < 40, and N > 40, respectively. The opacity, 
corrected for stimulated absorption, is then 

K * a = X H pN A a a A (l-J^)-\ (13) 

Since N n (N) = for N > 40, this absorption and emission process 
plays a role only during the very early phase of collapse. Typically, at 
densities near p ~ 10 12 g cm~ 3 , k* — > 0. Eq. 12 is only approximate and 
better estimates of the v e — A absorption cross section are in the offing 
(Langanke & Martinez-Pinedo 2003; Pruet & Fuller 2003). 
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3.4 Vi + p — > v>i + p: 

The total Vi—p elastic scattering cross section for all neutrino species 



is: 



a °( £u ^ ^4sin 4 ^-2sin 2 ^+ ) , (1-0 



p 4 \m e c 2 

where 9\y is the Weinberg angle and sin 2 9w — 0.23. In terms of C' v = 
1/2 + 2 sin 2 9w and C' A = 1/2 (note primes), eq. (14) becomes (Schinder 
1990): 

*> = T fe) 2 ^ " 1)2 + 39a{C ' a ~ 1)2 ] • (15) 

The differential cross section is: 
where 

(C^-l) 2 -^-!) 2 
dp (C(,-l) 2 + 3 5 2 (^-l) 2 • U0 

Note that S p , and <5„ below, are negative (5 P ~ —0.2 and 5 n ~ —0.1) 
and, hence, that these processes are backward-peaked. 

The transport (or momentum-transfer) cross section is simply 

< = y (^) 2 Wv - i) 2 + tf A (C A - if] . (is) 



where 



-f-/^(i-rt*- ff( (i-l« 



(19) 



3.5 i/j + n, — > i/i + n: 

The total v; L — n elastic scattering cross section is: 



The corresponding differential cross section is: 



^ = + PI) 



where 



(22) 
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The transport cross section is 

The fact that 5 P and 8 n are negative and, as a consequence, that cr' r is 
greater than Oi increases the neutrino-matter energy coupling rate for a 
given neutrino flux in the semi-transparent region. 

Horowitz (2002) derived expressions that include a weak magnetism/recoil 
correction analogous to those previously discussed for the charged-current 
absorption rates u e n <-> pe~ and v e n <-> ne + . We take the following form 
for the weak magnetism/recoil correction, a fit to the actual correction 
factor for the transport cross sections: 

<p < P (! + Cw M ^K lP ) , (24) 

where for neutrino-neutron scattering Cw M — —0.766, for neutrino- 
proton scattering Cjy M — —1.524, for anti-neutrino-neutron scatter- 
ing Cw M — —7.3656, and for anti-neutrino-proton scattering Cw M — 
-6.874. 

In fact, neutrino-nucleon scattering is slightly inelastic and when this 
is germane, as with and v T neutrinos, the more general formalism of 
§4, §5, §6 is necessary. 

3.6 Ui + A Vi + A: 

In the post-bounce phase, nuclei exist in the unshocked region exte- 
rior to the shock. At the high entropies in shocked protoneutron star 
atmospheres there are very few nuclei. There are alpha particles, but 
their fractional abundances are generally low, growing to interesting lev- 
els due to reassociation of free nucleons just interior to the shock only at 
late times. However, nuclei predominate on infall and neutrino-nucleus 
scattering (Freedman 1974) is the most important process during the 
lepton trapping phase. 

The differential Vi — A neutral-current scattering cross section may 
be expressed as: 

' h 1 (-^) 2 A 2 { W C FF + C LOS } 2 (S ion ) (1 + n) , (25) 



dQ 647r \m e c 2 
where 

97 

W = l-— (l-2sin 2 6V), (26) 

JTl. 

Z is the atomic number, A is the atomic weight, and (Si on ) is the ion- 
ion correlation function, determined mostly by the Coulomb interaction 
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between the nuclei during infall. (<Sj on ), in eq. (25) was investigated by 
Horowitz (1997) who approximated it with the expansion 



(Sion(e)) = 



1 + exp I 



where 



r 



(Ze) 2 1 
a kT 



hca 



1/3 



(27) 



(28) 



a is the interparticle spacing, rii on is the number density of ions, T is the 
ratio of the Coulomb potential between ions to the thermal energy in the 
medium, and are specified functions of T for each neutrino species. 

Leinson et al. (1988) have investigated the electron polarization cor- 
rection, Clos-, an d find that 



Clos = 



1 w 



1 + 4 sin 2 



1 + (kr D ) 2 



(29) 



where the Debye radius is 



7T 



h 2 c 



(30) 



k 2 = |p — p'| 2 = 2(£j,/c) 2 (l — ju), pf and Ep are the electron Fermi 
momentum and energy, and a is the fine-structure constant (~ 137 -1 ). 
Note that ro ~ Wh/pp in the ultra-relativistic limit (pf » m e c). The 
Clos term is important only for low neutrino energies, generally below 
~ 5 MeV. 

Following Tubbs & Schramm (1975) and ?, the form factor term, Cff, 
in eq. (25) can be approximated by: 



Cff 



,-y{l-n)/2 



(31) 



where 



H e »>^ysv) 2 (^) 



2/3 



and (r 2 ) 1 / 2 is the rms radius of the nucleus. Cff differs from 1 for large 
A and e u , when the de Broglie wavelength of the neutrino is smaller than 
the nuclear radius. 

When (Sion) = Cff = Clos + 1 = 1, we have simple coherent Freed- 
man scattering. The physics of the polarization, ion-ion correlation, and 
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form factor corrections to coherent scattering is interesting in its own 
right, but has little effect on supernovae (Bruenn & Mezzacappa 1997). 
The total and transport scattering cross sections for v j — a scattering 
(Z = 2;A = 4) are simply 

a Q = ^ = 4a (-^) 2 sin 4 ^. (32) 

4. The Neutrino Scattering Kernel and a 

Simplified Approach to Inelastic Scattering 

Many authors have studied inelastic neutrino-electron scattering as 
an important energy redistribution process which helps to thermalize 
neutrinos and increase their energetic coupling to matter in supernova 
explosions (Bruenn 1985; Mezzacappa & Bruenn 1993abc). Compara- 
tively little attention has been paid to inelastic neutrino-nucleon scatter- 
ing. Thompson, Burrows, & Horvath (2000) and Raffelt (2001) showed 
that, at least for and v T neutrinos, this process cannot be ignored. 
Here, we review the Legendre expansion formalism for approximating 
the angular dependence of the scattering kernel, detail our own imple- 
mentation of scattering terms in the Boltzmann equation, and include a 
discussion of neutrino-nucleon energy redistribution. In §5, we present 
an alternate approach involving dynamical structure factors that is more 
easily generalized to include many-body effects (§6). 

The general collision integral (RHS of the Boltzmann equation) for 
inelastic scattering may be written as 

CT u [fu] = (1 - /„) J ^0^f' u R hl (e u ,e' ul cos6) 

(33) 

= -scatt _ -scatty (34) 

where cos 6 is the cosine of the scattering angle, e u is the incident neu- 
trino energy, and e' v is the scattered neutrino energy. Although we sup- 
press it here, the incident and scattered neutrino phase space distribu- 
tion functions (f u and fl, respectively) have the following dependencies: 
fu = fu(r,t, fJ,,£ u ) and f' v = f' v {r, t, fj,', £' v ). [i and fx' are the cosines of 
the angular coordinate of the zenith angle in spherical symmetry and 
are related to cos 6 through 



cos 9 = nn'+ [(1 - n 2 )(l - n' 2 )] 1/2 cos(0 - (/)'). (35) 



Neutrino Opacities in Nuclear Matter 



11 



The only difference between f v and T v in §2 is that here f v has explicit 
fj, and e v dependencies. R m is the scattering kernel for scattering into 
the bin {e v , /j,) from any bin (e' u , ji') and R ont is the scattering kernel 
for scattering out of the bin (e u , fi) to any bin [e' v , //). The kernels are 
Green's functions that connect points in energy and momentum space. 
One may also write R(e u ,e' u ,cos6) as R(q,u>), where u>(= e v — e' u ) is 
the energy transfer and q(= [ef, + e'„ — 2e u e' u cos 9] 1 / 2 ) is the momentum 
transfer, so that the kernel explicitly reflects these dependencies (§5). 

An important simplification comes from detailed balance, a conse- 
quence of the fact that these scattering rates must drive the distribution 
to equilibrium. One obtains: R in = e _/3w i? out , where = 1/T. There- 
fore, we need deal only with R ont . The scattering kernels for inelastic 
neutrino-nucleon and neutrino-electron scattering depend in a compli- 
cated fashion on scattering angle. For this reason, one generally approx- 
imates the angular dependence of the scattering kernel with a truncated 
Legendre series (Bruenn 1985). We take 

°° 21 + 1 

i? out ( e; ,,4,cos0) = J]— ^$(e„,4,cos0)JKcos0), (36) 

1=0 

where 

r+l 

$ l (e u ,e' u ) = J d(cose)R° ut (e u ,e' l/ ,cose)P l (cos6). (37) 

Figure 1 shows R out (e u , e' u , cos 6) for ^-electron scattering as a function 
of cos 6 for various values of e' v at e v = 20 MeV and a representative 
thermodynamic point. In practice, one expands only to first order so 
that 

R ont {e u , 4, cos 6) ~ -<S> Q {e v ,e! v ) + -$i(e„, 4) cos 9. (38) 

Figure 2 depicts the various Legendre contributions for the same point. 
Figure 3 provides the corresponding final-state energy distribution func- 
tion for various values of e„, using only the / = term. Figures 4 and 
5 are similar to Figs. 1 and 3, but are for inelastic ^-nucleon scattering 
at a representative thermodynamic point. 

A simple approach to handling energy redistribution involves calcu- 
lating angle and energy moments and deriving source and sink terms 
(eqs. 33 and 34; Thompson, Burrows, & Pinto 2003). Substituting into 
the first term on the right-hand-side of eq. (33) (the source) gives 



/•oo j / 12 f+l i-2-k 



1 3 

-$ + cos' 

(39) 
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Substituting for cos 6 using eq. (35) and using the definitions 

r+l 



1 f +1 



dufu 



and 

we have that 



i 



-i 



Hu= 2J 1 d ^ v 



scatt = ^ _ j j 



4lT 



c(2vr/ic) 3 J 



(40) 



(41) 



(42) 



Integrating over jjl to get the source for the zeroth moment of the trans- 
port equation, 



- l +l dufi sc&tt = — 

2j_ L ^ C (27Thcf 



poo 

/ de'^'J 
Jo 



^0^(1 - J„) - hiH u H' u 

' (43) 
Similarly, we can write the sink term of the Boltzmann equation collision 

term (second term in eq. 33), employing the Legendre expansion 



-scatt 
Xu 



47T 



c(2iThc) 3 J 

The contribution to the zeroth moment equation is then 



(44) 



Combining these equations, we find that 
r+l 4vr 



i$ (i - J' u )l - hiH„K 

(45) 



de' e' 2 
c{27rhc)3 J v v 



^0 



J' V {1 ~ Ju)e 



-ft* - (l _ j' v ) J^j _ - $1 ^^(e-^ 



(46) 

One can see immediately that including another term in the Legendre ex- 
pansion (taking R out ~ (l/2)$ + (3/2)$i cos<9 + (5/2)$ 2 (l/2)(3cos 2 0- 
1)) necessitates including P v and P^, the second angular moment of the 
neutrino phase-space distribution function, in the source and sink terms. 
While easily doable, we advocate retaining only the linear term. 
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5. Inelastic Neutrino Interactions with 
Relativistic Nucleons and Leptons 

We now explore more sophisticated formalisms for handling inelastic 
scattering processes in nuclear matter. In this section, we address the 
non-interacting nucleon case. 

For neutrino energies of interest to supernova, which are less than 
a few hundred MeV, we may write the neutrino-matter interaction in 
terms of Fermi's effective Lagrangian 

v l + B 2 ^l + B i (47) 

v l + B 2 ^v l + B 4 , (48) 

where Gf — 1.436 x 10~ 49 erg cm~ 3 is the Fermi weak coupling constant. 
When the typical energy and momentum involved in the reaction are 
small (compared to the mass and 1 / size of the target particle) the lepton 
and baryon weak charged currents are: 

In = V ; «7m( 1 -75)Vv, jw = (gv - 9Alb) V>2 ■ (49) 
Similarly, the baryon and neutrino neutral currents are given by 

= i>uly. (1 - 75) , fz = ^47 M (cy - ca75) ^2 , (50) 

where 2 and 4 are the baryon (or electron) initial state and final state 
labels, respectively (these are identical for neutral-current reactions). 
The vector and axial- vector coupling constants (cy,<?y & ca,9a) are 
listed in Table 1 for the various charged- and neutral-current reactions of 
interest. The charged-current reactions are kinematically suppressed for 
v a and v T neutrinos. This is because their energy E Ufi / UT ~ T < m M , m T . 
On the other hand, neutral-current reactions are common to all neutrino 
species and the neutrino-baryon couplings are independent of neutrino 
flavor. Neutrino coupling to the lepton in the same family is modified 
since the scattering may proceed due to both W and Z exchange; the 
couplings shown in Table 1 reflect this fact. 

From the structure of the current-current Lagrangian, we can calcu- 
late the differential cross section for neutrino scattering and absorption. 
We are generally interested in calculating scattering/absorption rates in 
matter. Hence, it is convenient to express results in terms of the differ- 
ential scattering/absorption rate. For a neutrino with energy E\, this is 
given by dT(Ei) = J2i c da(Ei)i/V , where da is the differential cross 
section, the sum is over the target particle in volume V, and c = 1 is the 



pec 



L -int 



Gf , 

Gf .„ p 
v/2 » 3z 



for 
for 
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Reaction 


<V,Sv 


CA,gA 


F-2 


Vi + n^ Vi + n 

Vi + P — > Vi + p 

v e + e~ — ► v e + e~ 
v^/t + e~ — » Vy.iT + e~ 
i/ e + n <-» e~ + n 


-1/2 
1/2 - 2 sin 2 6W 
1/2 + 2 sin 2 W 
-1/2 + 2sin 2 (9w 
1 


-9a/2 
9a/2 
-1/2 
-1/2 

9A 


— 2fi n sin 2 9w — S 
— 2 fi p sin 2 9w + 5 





28 



Table 1. Neutral and charged current vector, axial and tensor coupling constants. 
The axial coupling qa = —1.23, the weak mixing angle sin 2 (9^=0.23. fi n = —1.913 
and Hp — 1.793 are the anomalous magnetic moments ol the neutron and proton, 
respectively and 8 — (/j, p — n„)/2 = 1.853. 



(relative) velocity of the neutrinos. For reasons that will soon become 
clear, it is useful to express the differential rate in terms of the current- 
current correlation function for the target particles (also called medium 
polarization tensor and is related to the dynamic structure factor) rather 
than the S-matrix element involving free asymptotic states. 

The differential scattering rate of matter is given by (Horowitz & 
Wehrberger 1991; Reddy et al. 1998): 



d 2 Q 3 dE 3 



£3 

32vr 2 E[ 



[1-/3(^3) 



exp ( zf^ztoiT ) 



Tm 



af3yrR 



L a/J n 



a3 



(51) 



where the incoming neutrino energy is Ei, the outgoing neutrino (or 
electron) energy is E3, and uj = E\ — £3. The factor [1 — exp((— to — [12 + 
H4)/T)] arises due to the fluctuation-dissipation theorem, since particles 
labeled '2' and '4' are in thermal equilibrium at temperature T and in 
chemical equilibrium with chemical potentials ^2 and ^4, respectively. 
For neutral-current processes, \i± = (12, and for charged-current pro- 
cesses, such as u e + n — > e~ + p, ^2 = and ^4 = fj, p . The final-state 
Pauli blocking of the outgoing lepton is accounted for by the factor 
(1 — fs(E 3 )). The lepton tensor L a p is given by 

L af3 = %[2k a k p + (k ■ q)g af3 - (k a q^ + q a kP) T ie a ^ y Wq u ] , (52) 



where the sign of the last term is positive for anti-neutrinos and negative 
for neutrinos. 

The medium or target particle retarded polarization tensor is given 

by 



Imng /3= tanh( " + ( ^-^ )lmn a ,, 



(53) 
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where II Q/ g is the time-ordered or causal polarization and is given by 

= "* / Tr P^Cp) r « G 4 (p + q) Fp)} . (54) 

Above, is the incoming neutrino four-momentum and is the four- 
momentum transfer. In writing the lepton tensor, we have neglected the 
electron mass term, since typical electron energies are of the order of a 
hundred MeV at high densities. The Green's functions Gi(p) (the index 
i labels particle species) describe the propagation of baryons at finite 
density and temperature. In a Fermi gas, the fermion Green's function 
is given by 

G ^ = luT ~ ' < 55 ) 

p + /U7o - m 

where \i is the chemical potential and m is the mass. For the neutral 
(charged) current, the vertex function is cv{gv)l^ f° r the vector 
current and ca(9a)7ix75 for the axial current. Given the structure of the 
particle currents, we have 

U af3 = c 2 v U v ^ + c\U% + 2c v c A U v J. (56) 

For the neutral-current, the vector polarization 11^ ^ is obtained by sub- 
stituting {r^r^} :: {7^,7,3}; the axial polarization is obtained 
by substituting {T a ,Tp} :: {7^75,7/375}; and the mixed part IT^| is 
obtained by substituting {T a ,Fp} :: {j a "f 5 ,lp}- 

Using vector current conservation and translational invariance, Il^g 
may be written in terms of two independent components. In a frame 
where q^ = (to, \q\,0, 0), we have 

Q 2 

n T = n& and n L = -^< . 

\Q\ 

The axial current-current correlation function can be written as a vector 
piece plus a correction term: 

n^n^ + ^ir 4 . (57) 

The mixed, axial-current/vector-current correlation function is 

IL™=ie m0 q a IL VA . (58) 
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These simplifications enable us to write the differential rate in eq. (51), 
as follows 



rfr £3 2 

x [A Sidg^u;) + S 2 (|«1,w) ±£S 3 (|<M], 
4EiS 3 + ql 



A 



2q 2 



B — Ei + £3 



(59) 



where the plus is for neutrinos and the minus is for anti-neutrinos. The 
three response functions, Si, S 2 , and S3 are given by 



S2(\q\,v) = 
S3(\q\,v) = 



(c v 2 + c A 2 ) (Imng( g ,^)+Imng(g,cc>)) 

1 - exp (-P(u + H2- A*4)) 
(cy 2 + c A 2 ) Im Ii$(q,u) - c\ Im n^(g,(j) 
1 - exp (-/?(w + /U 2 - ^4)) 
2cyCA Im n^ A (g,w) 



1 - exp (~P(uj + fjL 2 - A*4)) 



(60) 
(61) 
(62) 



The imaginary parts of the polarization functions needed to evaluate 
the response functions have been calculated explicitly by Horowitz &; 
Wehrberger (1991,1992) and Reddy et al. (1999). Neutrino scattering 
and absorption kinematics probes space-like kinematics where q 2 < 
(this is true only for massless neutrinos and electrons; however, when 
the typical lepton energies are large compared to m e , restriction to space- 
like kinematics is a good approximation). For space- like excitations, 



Im n£(g,w) 

Im U§(q,uj) 

Im IL%(q,u) 
Im n$ A (q,uj) 



2ir\q\ 
q 2 

4-7T 1 I ' 

,2 



2vr|g| 
5^^0 + 2/1], 



(63) 

(64) 

(65) 
(66) 



where Mj is the rest mass of particle i. The one-dimensional integrals 
In = tanh ^ + (^~ fH^ J°° dE En [F{E ^ E + UJ)+F{E + W) E)] ; 

(67) 
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where F(x,y) = /2(x)(l — f/±{y)) and fi(E) are Fermi-Dirac particle 
distribution functions. The lower limit on the integral, e_, arises due to 
kinematical restrictions and is given by 



where 

q A — uj z 

These integrals may be explicitly expressed in terms of the Polylogarith- 
mic functions 

Li n {z)= Lln -^ dx, Lii(s) = ln(l-x). (70) 
Jo x 

This Poly logarithm representation is particularly useful and compact: 

Jo = T z (l + ^ V (71) 

(73) 

where z = (u + (^2 — ^a))/T and the factors £ n are given by 

in = Li n (-ai) - Li n (-a 2 ) , (74) 

with 

«i = exp ((e_ - H2)/T) , 0:2 = exp ((e_ + a; - Ha)/T) . (75) 
In the case of neutral currents, some of the terms above simplify 



z = -, 1*2 = 1*4, e_ = -- + -Wl -4— |- . (-()) 



The total scattering rate is the double integral in (q,u) space: 

p2 i-Ei r2Ei—u 

= dw(l-f 3 (E 3 )) dqqq 2 n [AS 1 +S 2 ±BS 3 ]. 



(77) 
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Eq.(77) allows us to calculate the cross section per unit volume, or equiv- 
alently the inverse mean free path. This naturally incorporates the ef- 
fects of Pauli blocking, Fermi and thermal motion of target particles, and 
the contribution of relativistic terms to the electron and baryon currents. 
The formalism reviewed in this section provides a unified description of 
reactions (both charged and neutral currents) involving both nucleons 
and electrons. In the absence of strong and electromagnetic interactions, 
this is almost all there is to the subject of neutrino opacities in matter 
with nucleons and leptons. The only omission here is the effect of weak 
magnetism (§3.1). This arises due to tensor coupling between neutrinos 
and nucleons and is discussed next. 

The vector and axial-vector structure of the weak hadronic currents 
given in eqs. (49) and (50) are valid only at low energy {E v <C M). 
When the neutrino energy becomes comparable to the intrinsic energy 
scale that characterizes the hadron (its size and mass) additional con- 
tributions arise. These include form factors and the tensor (weak mag- 
netism) couplings. To leading order in the four-momentum transfer ex- 
pansion (q^), the hadronic neutral and charged weak currents take the 
form 

& = ^A[gv{ql)^ + iF 2 {ql)^^-g A {ql)^ lb ]^ (78) 

fz = ^ [M^-iFiiq 2 ^^- c\(ql)n,] fh- (79) 

The q 2 dependence of the vector and axial couplings is weak and one 
can safely ignore these corrections, but the contribution from weak mag- 
netism is not negligible (Vogel 1984). The tensor coupling constants for 
the various reactions of interest are given in Table 1. For a neutrino 
of energy E u ~ 50 MeV, a typical \q\ ~ 2E V ~ 100 MeV, the weak 
magnetism correction to the vector current is then F2\q\/(2M) ~ 6%. 
Perhaps, more importantly, these corrections appear with opposite sign 
for neutrinos and anti- neutrinos (Vogel 1984; Horowitz 2002), resulting 
in 10 — 20% differences in the relative cross sections. 

Tensor coupling introduces new structure to the current-current cor- 
relation function given in eq. (56) and corresponding changes to the 
response functions given in eq. (62). These corrections have been com- 
puted in recent work by Horowitz & Perez-Garcia (2003). To order 
q 2 /M, the tensor- vector, tensor-axial vector and tensor-tensor polar- 
ization functions contribute to the neutrino response. These are ob- 
tained from eq. (54), with the following substitutions: for the tensor- 
vector, T Q = iF2(T a ,5q S /2M and Tp = c v jp; for the tensor-axial, T a = 
iF2<T a7 sq S /2M andT/3 = Cajpjs; and for the tensor-tensor, T Q = iF2(T a ^q 5 /2M 
and = iF2<Ji3 } sq S /2M. These polarization functions can be related to 
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the vector, axial, and vector-axial polarization tensors discussed earlier. 
Explicit forms for the imaginary parts of these polarization tensors can 
be found in Horowitz & Perez-Garcia (2003). 

6. Strong and electromagnetic correlations 
between nuclei, nucleons, and electrons 

The results of the previous sections provide a complete description of 
neutrino cross sections in an ideal Fermi gas of nucleons and leptons. 
Here, we address non-ideal corrections to neutrino scattering rates aris- 
ing due to interactions between nucleons, electrons and nuclei. The role 
of correlations in neutrino opacities was first appreciated in the pioneer- 
ing works of Sawyer (1975) and Iwamoto & Pethick (1982). They showed 
that strong interactions between nucleons could significantly modify (by 
as much as a factor of 2-3) neutrino cross sections at densities and tem- 
peratures of relevance to supernovae. In the past decade, there have 
been several attempts to incorporate these effects in calculations of the 
neutrino opacity of dense matter. While much progress has been made in 
understanding qualitative aspects of these corrections, quantitative pre- 
dictions for the neutrino cross sections at high density remain elusive. 
We present a pedagogic review of the qualitative findings of model calcu- 
lations, provide an overview of the current state of the art of many-body 
effects in neutrino cross sections, and comment on their shortcomings. 

6.1 Plasma of Heavy Ions: Neutrino Opacity at 
p ~ 10 12 g/cm 3 : 

The simplest system, with non-trivial many-body dynamics, which is 
relevant in the supernova context is a plasma of heavy nuclei (like Fe) 
immersed in a background of degenerate electrons. Freedman (1974) 
showed that the dominant source of opacity for neutrinos in such a 
plasma is the coherent scattering off nuclei. Low-energy neutrinos can 
couple coherently, via the vector neutral current, to the total weak charge 
of the nucleus. The weak charge of a nucleus with A nucleons and Z 
protons is given by Qw = (2Z — A — 4Z sin 2 6w)/2- 

For neutrino-nucleus scattering in a plasma, we must properly ac- 
count for the presence of other nuclei, since scattering from these dif- 
ferent sources can interfere. In the language of many-body theory, this 
screening is encoded in the density-density correlation function. It is 
therefore natural to express the scattering rate in terms of these cor- 
relation functions as described earlier in §5. To motivate the relation 
between the cross section and the density-density correlation function 
we begin by noting that the effective Lagrangian describing the neutral- 
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current interaction of low-energy neutrinos with nuclei is given by 

Gv 

L NC = Qw ^ f (80) 

where 1^ = V^^(l — 75)^ is the neutrino neutral current. Nuclei are 
heavy, and, correspondingly, their thermal velocities are small (v = 
yjT/M <C c). For simplicity, we assume that nuclei are bosons char- 
acterized only by their charge and baryon number. In this case it is 
an excellent approximation to write the neutral current carried by the 
nuclei as = ijpij) 5q. We can write the differential scattering rate in 
terms of the density operator in momentum space given by 

p(q,t)=^^= exp(iq-fi(t)), (81) 

i=l-N 

where the sum is over ./V particles in a volume V which are labeled 
i = 1 • ■ ■ N. The rate for scattering of a neutrino with energy E v to a 
state with energy E' v = E v — 00, at an angle 6, with a momentum transfer 
q is given by 

__ = _Z Q 2, (1 + cos9)£ - m „ ) . (82) 

The function S(\q\,tu) is called the dynamic structure function. It em- 
bodies all spatial and temporal correlations between target particles aris- 
ing from strong or electromagnetic interactions and is given by 

/oo 
dt exp(iwt) (p(q,t)p(-q,0)) . (83) 
-00 

In the above, (• • • ) denotes the ensemble average per unit volume. Note 
that S*(|g|,u;) is normalized such that in the elastic limit for a non- 
degenerate, non-interacting gas S(\q\,uj) = 2ir S(lo) n, where n = N/V 
is the density of particles. 

To calculate the structure function we need to solve for the dynamics 
(r(t)) of the ions as they move in each others presence, interacting via the 
two-body ion-ion interaction potential. The electrons are relativistic and 
degenerate, and the electron Fermi momentum kp = (3tt 2 Z n) 1 / 3 . The 
timescale associated with changes in their density distribution is rapid 
compared to the slow changes we expect in the density field of the heavy 
ions. Electrons almost "instantaneously" follow the ions and screen the 
Coulomb potential between ions. Correspondingly, there is small excess 
in the electron distribution around each ion. As first discussed by Lein- 
son et al. (1988), this excess will also screen the weak charge of the nu- 
cleus to which the neutrinos couple. Electronic screening was discussed 
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earlier in §3.6 and eq. (25). Here, we focus on the dynamics of ions inter- 
acting throug h the sc reened potential V(r) = Z 2 e 2 exp(— r / X e ) / (Airr) , 
where A e = \J 4a cm /7r kp 1 is the electron Debye screening length. 

At 10 12 g/cm 3 , the typical inter- ion distance is d <J 30 fm, A e > 80 fm, 
and the typical ion-ion interaction energy, E pot ~ Z 2 e 2 /(And), is large. 
For temperature in the range 1 — 5 MeV, the ratio of the potential energy 
to the kinetic energy, T, equals Z 2 a em /(d kT). For nickel-like nuclei with 
Z = 28, r = 8.8 for T = 4 MeV. The dynamics of such a strongly cou- 
pled plasma is not amenable to analytic methods of perturbation theory 
and approximate non-perturbative methods have had limited success in 
describing these systems. For a classical system of point particles inter- 
acting via a 2-body potential it is possible to numerically simulate the 
real-time dynamics of the system. Such simulations confine N particles 
to a box with periodic boundary conditions, calculate the force on each 
particle at any time, and evolve the particles by using their equations of 
motion (r = Force/M, where M is the mass of the ion). This technique, 
which goes by the name molecular dynamics (MD), has been used ex- 
tensively in condensed matter physics, plasma physics, and chemistry. 
For an early application of this technique to the study of the response 
of a one-component plasma see Hansen, McDonald, & Pollock (1975). 

The de Broglie wavelength of the ions is Ad — 4.6 x A" 1 f m > 
where A is the mass number of the ion and XM e v is the temperature 
in MeV. For T > 1 MeV, Ajy <C d, the ion-plasma is classical, and we 
may implement MD to calculate S(\q\,uj). Fig. (6) shows the results 
of such a calculation (from Luu et al. 2004). We chose to simulate 
54 ions (A=56,Z=28) in a box of length L = 200 fm. This corre- 
sponds to an ion density n = 6.75 x 10~ 6 fm~ 3 or the mass density 
p ~ 6 x 10 12 g/cm 3 . The background electron density was chosen to 
make the system electrically neutral. For a classical simulation, the sys- 
tem is characterized by the plasma V defined earlier. Figure (6) shows a 
comparison between the results for S(q,co) obtained by MD simulations 
(dots with error bars), the free (Boltzmann) gas response (dashed- line), 
and the response obtained using Random Phase Approximation (RPA). 
The response in RPA is discussed in the next section. The results are 
shown as a function of ^/^piasmon) where to is the energy transfer, and 
^piasmon = V ' ^ «cm n/M is the plasma frequency of the classical heavy- 
ion plasma. 

6.1.1 Random Phase Approximation (RPA). We motivate 
and present a heuristic derivation of the structure function calculated 
using the random phase approximation. To begin, we note that the dy- 
namic structure factor for a non-relativistic and non-interacting Maxwell- 
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Boltzmann gas with a chemical potential ha and number density n is 
given by 



Smb (q, oj) = J ^3 f(E p ) 2vr S (E p - - u) 
n /27T M . ^ . 

= ^ V "It" exp( - e +> ■ < 84 > 



This follows directly from the equation of motion of the free particles 
and eq. (83). When an external current couples to the interacting sys- 
tem at wavelengths large compared to the inter-particle spacing, scat- 
tering amplitudes from different particles can interfere. In the case of 
neutrino scattering this leads to screening of the weak charge. RPA is 
known to provide an adequate description of this dynamic screening in 
the long-wavelength limit. Pioneering work by Bohm & Pines (1982) 
established RPA as a useful non-perturbative approximation. It is used 
widely in nuclear and solid state physics to describe the long-wavelength 
response of many-particle systems. In particular, RPA accounts for dy- 
namic screening and is able to correctly predict the presence of collec- 
tive modes. These correspond to a plasmon in the case of unscreened 
Coulomb interactions and phonons in the case of screened interactions 
(as in the present case). The RPA results in Fig. 6 clearly show the 
presence of these collective states. 

Screening in RPA is accounted for by the selective re-summation of 
an infinite series of bubble graphs. The bubbles are free excitations and 
are connected by the two-body interaction potential. The lower diagram 
in Fig. 7 represents this series. The solid lines represent propagation of 
ions, the wavy line represents the external current (the neutral current), 
and the dashed lines represent the screened Coulomb interaction. The 
diagram illustrates screening of the coupling between the external probe 
(neutrinos in the present case) and particles in the plasma through the 
intermediate particle-hole excitations. The bubble-sum is a geometric 
series in the interaction and can be evaluated analytically for simple 
interactions. For the ion-plasma, the dynamic structure function in RPA 
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is given by 



<Srpa(<7» = z 2 f „ z r Im 

1 — exp (—p tvjj) 



no(g,o;) 
1- V c (g) n (g» 



,(87) 



where F c (g) = a cm Z 2 ( ^l^ ) ■ ( 88 ) 

The expression for Vc(q) is valid for q <C kp, where A e is the Debye 
screening length and k-p is the Fermi momentum of the background de- 
generate electron gas. The free-gas polarization function 

IW, w) = i f ^ G(p) G(p + q) , (89) 

where G(p) = (po — (p 2 /2M — ha)) -1 is the free particle, non-relativistic, 
Green's function and fiA is chemical potential of the ions. Analytic 
expressions for the real and imaginary parts of Ilo were first obtained by 
Lindhard (1954). In the case of a Boltzmann gas these are given by 

6 M 2 

Imn (g» = tL (l-exp(-/?M) exp(-Q 2 ,) , (90) 

Z7T q 

t> rr r- \ t> f duJ> Imn o(g,w') . . 

Ren (g,w = V / — , 91 

J 7T LO — OJ' 

where Q\ was defined earlier in eq. (85). 

The comparison in Fig. 6 shows that both RPA and the free-gas 
response differ quantitatively. MD is exact in the classical limit. Cor- 
rections to classical evolution are governed by the smallness of the ex- 
pansion parameter Ap/d. In the present case (A^/d) <S 1%. In the 
left panel, the response for long wavelengths (\q\ = 2tt/L ~ 6 MeV) is 
shown. In this case, the discrepancy between RPA and MD is not signifi- 
cant, supporting the expectation that RPA provides a fair description at 
long wavelengths. It captures, albeit much too sharply, the phonon peak 
(which in solid state parlance is called the Bohm-Staver sound) seen in 
MD simulations. The right panel shows results for q = 6ir/L ~ 18 MeV. 
This corresponds to a wavelength that is comparable to the inter-ion 
distance. Here, the plasmon mode is damped in RPA by single-pa rticle 

excitations, since tJpl asmon ~ (/^thermal, where ^thermal = t/3 kT/M is the 

thermal velocity of the ions. Nonetheless, this damping is still too weak 
to dissolve the plasmon peak. In contrast, results obtained using MD do 
not show any collective behavior, suggesting that RPA underestimates 
damping. 

These results clearly illustrate two important features of the response 
of strongly-correlated systems: (1) correlations can greatly affect the 
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shape of the low-energy response; and (2) approximate many-body meth- 
ods such as RPA can capture qualitative aspects of the response at large 
wavelength, but do poorly even at moderate wavelength. The failure 
of RPA can be attributed to the neglect of collisional damping from 
the excitation of multi-particle states. As the ions move, they "simul- 
taneously" scatter off several correlated neighboring particles. In RPA, 
these correlations are neglected. Ignoring these multi-pair excitations is 
a reasonable approximation in the long- wavelength limit. We will jus- 
tify this in our subsequent discussion of sum rules (§6.1.2), where we will 
argue that the spectrum of multi-pair excitations vanishes as q 2 in the 
q — > limit. However, the results indicate significant contribution from 
multi-pair excitations even at moderate wavelength. Finding ways to 
incorporate collisional damping of single-particle motion in the response 
has a long history in solid state physics, starting with the pioneering 
work of Kadanoff & Martin (1963). We will not review these techniques 
here. Our focus was to illustrate the importance of many-particle dy- 
namics and provide a baseline for comparison. The lessons learned here 
will prove valuable in gauging the validity of results that will discussed 
in §6.3, where we evaluate the results obtained using RPA. 

6.1.2 Sum Rules. The dynamic structure factor is subject to 
constraints arising due to thermodynamic considerations and symmetries 
of the Hamiltonian (conservation laws) describing the many-particle sys- 
tem. These constraints, called sum rules, provide valuable guidance in 
developing approximation schemes needed to compute S(\q\,uj). In the 
following, we discuss two such sum rules: (i) the compressibility sum 
rule; and (ii) the energy-weighted sum rule. 

The compressibility sum rule is related to the long-wavelength limit 
of the static structure function which is defined as 

1 POO f J . 

In the long-wavlength limit (q — ► 0), the static structure function is 
related to the isothermal compressibility by the following relation 

lim Si* = — — lim (p(-q,0)p(q,0)) 

= nkT IC T , (93) 
1 fdV\ 

wh ere K T = -- (-^ (94) 

is the isothermal compressibility and n is the particle number density. 
The first of the above relations follows directly from the definition of 
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S(q,co) in terms of the density-density correlation function in eq. (83). 
The above relation permits us to verify consistency between the neutrino 
response functions and the equation of state. 

The energy-weighted or F sum rule for any operator O is given by 

[ ( izfghgM ) s (m = ( [w,6m,e>w] > (95) 

In the absence of any velocity-dependent terms in the Hamiltonian, and 
for the specific case of the density operator, the right hand of eq. (95) 
is independent of interactions and is given by 

( [[H,p(q)],p(q)} )=n£j. (96) 

Note that the F sum rule vanishes in the long-wavelength limit when the 
operator p(0) commutes with the Hamiltonian. This is a consequence 
of the conservation of particle number and provides a useful constraint 
for the large u behavior of S(q,co) when it is computed using approxi- 
mate many-body methods. For the response of the ion-plasma discussed 
above, the results obtained in the case of the free Boltzmann response, 
the RPA, and molecular dynamics satisfy this F sum rule. Despite the 
large differences in the low-energy strength. Another important con- 
sequence of the F sum-rule is that it requires response at long wave- 
lengths to vanish as q 2 . This implies that multi-particle excitations are 
suppressed at long wavelength when the operator commutes with the 
Hamiltonian. It also offers an explanation for the trend seen in Fig. 6, 
where RPA results improve at large wavelength. 

6.2 Pasta Phases: Neutrino Opacity at p > 10 13 
g/cm 3 

With increasing density, the nuclei get bigger and the internuclear 
distance becomes smaller. Under these conditions, the nuclear surface 
and Coulomb contributions to the free energy of the system become 
important. Ravenhall, Pethick, & Wilson (1983) showed that it is ener- 
getically favorable for nuclei to deform and assume novel, non-spherical, 
shapes such as rods and slabs. Further, the energy differences between 
these various shapes are small AE ~ 10 — 100 keV. The dynamics of 
such an exotic heterogeneous phase is a complex problem involving sev- 
eral energy scales. Both strong and electromagnetic interactions play 
an important role. For temperatures of interest, T <J 5 MeV, the de 
Broglie wavelength and the inter-particle distance are comparable and 
quantum effects cannot be neglected. Recently, Watanabe et al. (2003) 
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have studied the behavior of such matter using the techniques of quan- 
tum molecular dynamics and also find rod- and slab-like configurations. 
How does the heterogeneity and existence of several low-energy excita- 
tions involving shape fluctuations influence the response of this phase 
to neutrinos? In the simplest description, the structure size (r) and 
the inter-structure distance (R) characterize the system. We can expect 
that neutrinos with wavelength large compared to the structure size, but 
small compared to the inter-structure distance, can couple coherently to 
the total weak charge (excess) of the structure, much like the coherence 
we discussed in the previous section. The effects of this coherent en- 
hancement in the neutrino cross sections has recently been investigated 
by Horowitz, Perez-Garcia, & Piekarewicz (2004). In agreement with 
our naive expectation, their study finds that the neutrino cross sections 
are greatly enhanced, by as much as an order of magnitude, for neutrinos 
with energy 1/r > E v > 1/R. 

6.3 Nuclear liquid: Neutrino opacity at p > 10 14 
g/cm 3 

When the density of matter exceed 10 g/cm 3 , matter is expected to 
form a homogenous liquid of nucleons and leptons. To account for the 
effects of strong and electromagnetic correlations between target neu- 
trons, protons and electrons we must find ways to improve H a /3, the 
polarization tensor of free Fermi gas described earlier in eq. (54). This 
involves improving the Green's functions for the particles and the asso- 
ciated vertex corrections that modify the current operators. In strongly 
coupled systems, these improvements are notoriously difficult and no 
exact analytic methods exist. With few exceptions, most investigations 
of many-body effects in the neutrino opacities have employed the mean- 
field theory to improve the Green's functions. 

Dressing the single particle Green's functions must be accompanied 
by corresponding corrections to the weak interaction vertex function. 
The Green's functions in the mean field (Hartree) approximation in- 
volves the Schwinger-Dyson summation of the diagrams shown in Fig. 
7. The dressing of Green's functions implies the presence of an inter- 
action "cloud" surrounding the quasi-particle. This cloud modifies the 
coupling of an external probe to the quasi-particle. In §6.1, we discussed 
this screening of the weak charge due to Coulomb interactions between 
ions within RPA . Here, we consider screening due to both and strong 
and electromagnetic interactions. The modification of the coupling of 
an external current to the quasi-particles is shown in the lower panel of 
Fig. 7, where the solid squares represent interactions between nucleons. 
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6.3.1 Non-Relativistic RPA. At sub-nuclear density, nucleons 
are non-relativistic. In the non-relativistic limit the structure of the 
hadronic current greatly simplifies and the vector and axial currents can 
be directly related to the density and spin-density operators. Neutrino 
scattering off a non-relativistic neutron liquid was first investigated by 
Iwamoto & Pethick (1982). In this case, the scattering rate can be 
written as 



d cos 9dE' u 4vr 2 



x (c v 2 (l + cos#) S(\q\,to) + ca 2 (3-cos#) S A (\q\,u)) , 

(97) 

where the dynamic structure function for density fluctuations is 

/oo 
dt exp(iut) ( p(q, t)p(-q, 0) ) , (98) 
-oo 

and the dynamic structure function for spin-density fluctuations is 

/oo 
dt exp(iut) Stj ( CTi(q,t)(j(-q,0) > . (99) 
-oo 

The thermal ensemble average is denoted by (• • • ). It samples the opera- 
tors over all states of the many-particle system with weight exp(—7i/kT), 
where 7i is the Hamiltonian describing the neutron liquid. The formal- 
ism is easily extended to include protons (Sawyer 1975; Horowitz & 
Wehrberger 1991; Burrows & Sawyer 1998,1999). In this case, the dif- 
ferential neutrino-nucleon scattering rate from a liquid of neutrons and 
protons is given by, 

{Ait 2 )- l G 2 F E 2 [l - U(E 2 )} \(1 + cose)W) 2 S nn (q,uj) 



dw d cos 6 

+ (3 - cos#) 5 i [S^{q,u) + S£ n (q,u) - 2S£ n (q, u)] , (100) 

where Ei=E\ — uj. 

Following Burrows & Sawyer (1998), we consider a simple nuclear 
potential, with only central s-wave interactions, given by 

V = Vi + V2 T\ ■ r 2 + "^3 01 • <T2 + V4 01 • 2 ti • r 2 , (101) 

where a and r are 2x2 Pauli matrices acting in spin and iso-spin space, 
respectively. The structure functions, S (Fermi) and S A (Gamow- Teller; 
axial), are elements of separate 2x2 symmetric matrices. For the vector 
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dynamic structure function, S, we have 

\bp n (q,u) b nn (q,u).J 
The structure function matrix is given by, 



S(q,u) 



21m 



n(°)(g,w)[l - viq^iq,^)}- 1 ] (1 - e _/3w 



(102) 



where 



n(°)( g , w ) 



nf(g,o;) 




ni° ) (g,o;) / 



and nj, 0) and ni° } are given by the Fermi gas polarization functions and 
evaluated with the proton and neutron chemical potentials, respectively. 
The non-relativistic polarization for a Fermi gas is given by 



lmn (0) {q,uj) 
ReU {0 \q,uj) 



m 



2irqf3 



log 



l + e 



1 + e 



-Ql+0n-0w 



/die 
7T 



dJ Imn(°)(g, 



uj — ur 



where 



Q± = 

The potential matrix is, 



q 2m J 



(103) 
(104) 

(105) 



,2 \-l 



vi+v 2 + 4ire 2 (q 2 + q% F ) 
vi - v 2 



Vl + V 2 



where vi = J d?r Vi, and the term containing qxF is the Thomas- Fermi 
screened Coulomb potential {q^F = 4e 2 7r 1 / 3 (3n p ) 2 / 3 ). In Fermi liquid 
theory (FLT), the underlying interaction between (quasi)particle-hole 
states is directly related to thermodynamic quantities such as the sus- 
ceptibilities and specific heat (see Baym & Pethick (1991) for a pedagogic 
introduction to FLT). These effective interactions can differ greatly from 
the free-space nucleon-nulceon potentials. In the absence of experimen- 
tal determinations of the thermodynamic susceptibilities, microscopic 
theories of nuclear interactions are need to obtain these effective inter- 
actions. For the simple case of only central interactions, the suscepti- 
bilities of nuclear matter are encoded in four constants called Landau 
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parameters (Fo, Fq, Go, G' correpond to baryon density, isospin density, 
spin-density and spin- isospin density fluctuations, respectively). Fol- 
lowing Burrows & Sawyer (1998,1999) and for simplicity, we use these 
Landau parameters to determine the potential matrix, in lieu of a more 
developed nuclear interaction model. The parameters of the potential 
matrix are related to the Landau parameters by the following relations 

F +F' F -F' 
Vl+V 2 = -j^ V!-V2 = 

■ Go+G n Go— G n ' 

where No = 2m* kp/^ 2 is the density of states at the Fermi surface of 
nuclear matter, m* is the effective nucleon mass, and fcp is the nucleon 
Fermi momentum. Taking m* = 0.75m n as our fiducial value for the 
effective mass, we use parameters from Backman, Brown, & Niskanen 
(1985) and Brown & Rho (1981): F = —0.28; Fq = 0.95; G = 0;G' = 
1.7, obtaining, 

Vl = -7.4 x 10~ 6 MeV~ 2 
v 2 = 2.5 x 10" 5 MeV~ 2 
v 3 = 

v A = 4.5 x 10~ 5 MeV~ 2 . (106) 

For other values of the effective mass, we keep these potentials at the 
same value, which is to say we assume that the Landau parameters are 
proportional to m*/m. 

The form for the Gamow-Teller matrix, S A (q,co), is the same as that 
for S, except that the potential matrix is replaced by v A 

A = A>3 + V4 V 3 -V4 
~ \V 3 -V4 V 3 + V4 

Taking the matrix inverses leads to the following forms for the com- 
binations of structure functions that appear in eq. (100) 

S nn (q,cj) = 2Im[nflV] (1 - e"^)" 1 , (107) 

where 

D v = l-(v 1 + v 2 )U^ - ( Vl - v 2 ) 2 U^U^Q v - 1 . (108) 
Qv is given by the expression: 

Qv = l~ W(g 2 + (Ztf)- 1 !!^ - ( Vl + v 2 )U^ . (109) 
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If, as in eq. (106), we take = 0, we obtain the simple result for the 
axial-current terms, 



S A (q,co) = 21m 



(0), 



(1 



(110) 



For the Fermi term, since Cy = 1/2 — 2 sin 2 6w ~ 0, we drop the pro- 
ton structure function in eq. (100). Furthermore, we use the potential 
parameters given in eq. (106), and in eq. (108) we drop the third term. 
This term would have been significant had it not been for the Coulomb 
term in the denominator, an illustration of the importance of the explicit 
inclusion of Coulomb forces, even for the neutron density correlations. 
Since the i>jS are all real, we obtain for the structure factors used in eq. 
(100), 



where 
and 



C v = (1 - v F ReU^) 2 + 4(Im4°)) 2 , 



(111) 
(112) 



S A {q,u>) = 2[lml4°)(g, W ) + lmU^(q,u)} (1 - e^)" 1 ^ 1 , (113) 



where 

Ca = Cai + Ca2 ■ 
Cai and Ca2 are given by the expressions: 



Cm = [l - v GT (ReIlW(q,uj) + ReU^(q, 



and 



C A 2 = v 2 GT [lml4°) (q, uj) + Imlif (q, u) 



(114) 

(115) 
(116) 



The F in Sp(q, to) and the A in Sa(q, w) stand for Fermi and Gamow- 
Teller (axial) and vf and vqt equal (v\ + V2) and U4, respectively, in 
Fermi Liquid Theory. Sa{<1,u) in eq. (113) is now the entire axial term 
in eq. (100). C\>a is the correction factor due to many-body effects for 
a given momentum transfer (or scattering angle) and energy transfer. A 
similar procedure is employed for calculating the many-body corrections 
to the charged-current rates (Burrows & Sawyer 1999). Figure 8 portrays 
Sa(q,w) versus uj/q for representative parameters and clearly indicates 
the resonances and Cerenkov kinematics. 
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6.3.2 Multi-pair excitations and the Tensor Force. The 

simple model, based on Fermi Liquid Theory and the Random Phase Ap- 
proximation with central interactions, is rudimentary and relies on the 
poorly constrained determinations of the Landau parameters. While we 
can expect this simple analysis to capture qualitative aspects of many- 
body correlations it is not suitable for quantitative predictions. It has 
two important shortcomings, namely the neglect of multi-pair excitations 
and non-central interactions such as the tensor force (pion exchange), 
which is known to be important in nuclear systems. The discussions 
in §6.1 regarding the response of the ion-plasma clearly emphasized the 
need to incorporate finite propagation lifetimes for quasi-particles due 
to collisional damping. The role of collisional damping in the low-energy 
response is particularly well studied for the case of low-energy photon 
production from bremsstrahlung in many-particle systems and is called 
the Landau- Pomeranchuk-Migdal effect. Similar effects were first inves- 
tigated in the context of the response of nuclear matter by Raffelt & 
Seckel (1995). They found that the rate of spin fluctuations (l/r sp i n ) 
due to nucleon-nucleon collisions in the medium is rapid compared to 
the typical energy transfer to in neutrino scattering. Consequently, the 
axial charge is dynamically screened for small uj <J T spui , resulting in a 
suppression of the low-energy axial response. 

The redistribution of response strength in energy is a generic feature of 
many-particle systems arising due to the finite lifetime of quasi-particles. 
However, the situation in nuclear systems is unique due to the presence 
of a strong tensor force. This has recently been clarified by Olsson & 
Pethick (2002). The evolution of nucleon spin is dominated by tensor 
interactions, especially because the nucleon spin operator a does not 
commute with the tensor operator in the nuclear Hamiltonian. The F- 
sum rule, discussed earlier in §6.1, for the spin response function does 
not vanish in the long- wavelength limit, since [H tensor , o - ] / 0. Fur- 
ther, the relationship between the spin susceptibility and the Landau 
parameters is modified due to presence of the tensor interaction (see 
Olsson & Pethick (2002) for these revised relations). The preceding dis- 
cussion indicates that interactions that do not commute with the spin 
operator may be especially important in determining the low-energy re- 
sponse. Olsson & Pethick (2002) estimate that as much as 60% of the 
low-energy axial response at long wavelength may reside in multi-particle 
excitations. This preliminary estimate warrants further investigation re- 
quiring both the inclusion of the tensor force and multi-pair excitations 
in the axial response. 

When multi-pair excitations become important it is appropriate to 
work in terms of a correlated basis states rather than single particle 
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bare neutron and proton states. The correlated basis states are expected 
to be close to the energy eigenstates of the system. Consequently, in 
this basis the residual interactions are weak. Further, the ground and 
excited states in the correlated basis states contain multi-particle hole 
states and the quasi-particles of the correlated basis are superpositions 
of neutron and proton states of both spins. This is particularly relevant 
for nuclear matter, where pion exchange can transform both spin and 
isospin of the bare nucleons. Recently, Cowell & Pandharipande (2002) 
have computed weak interaction matrix elements in the correlated basis 
obtained using a two-body cluster expansion. They find that spin and 
iso-spin correlations play an important role and result in quenching the 
weak interaction transition rates by 20 — 25% at low energy. 

6.3.3 Relativistic RPA. The RPA polarization tensors that en- 
ter the neutrino scattering and absorption rates in a relativistic frame- 
work have been computed previously by Horowitz & Wehrberger (1992) 
and Reddy et al. (1999). The specific form of the polarization tensor in 
relativistic random phase approximation (RRPA) follows from the struc- 
ture of the bubble sum shown in Fig. 7. For a one-component Fermi 
system, such as neutron matter, the polarization functions appearing in 
eq. (51) are replaced by 

ttRPA ttMF i ttRPA n ttMF 

where D a p, is the relativistic equivalent of the potential matrix defined 
in §6.3.1. It describes the interaction between neutrons and IImf is the 
polarization tensor in eq. (66), but with the Green's functions for neu- 
trons computed in the mean-field approximation. In relativistic mean- 
field theories (RMF), which are inspired by the Walecka model, nucle- 
ons interact via the exchange of (fictitious) scalar, isoscalar-vector, and 
isovector- vector mesons (for recent review of RMF models see Serot &: 
Walecka 1997). The structure of D a p and the mean-field Green's func- 
tions in RMF theory may be found in Horowitz & Wehrberger (1992). 

The neutrino mean free paths computed in the context of these models 
indicate a moderate suppression of 10-20% in the neutrino cross sections 
at nuclear density. Despite the strong coupling between nucleons, the 
suppression is modest. This is because the interactions in RMF mod- 
els primarily affect the vector response, while the neutrino scattering is 
dominated by the axial response. The RMF model is constrained by its 
mean field predictions for the empirical properties of spin saturated mat- 
ter. It is necessary to supplement the RMF model with additional spin 
dependent interactions such as a pion exchange which do not contribute 
to the ground state properties in the Hartree or mean field approxima- 
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tion. However, these interactions modify the RRPA response. Model 
calculations in the RRPA which include spin-isospin dependent inter- 
actions (pion exchange and phenomenological short-range interactions) 
have been performed by Reddy et al. (1999). The results of these studies 
indicate that neutrino mean free path computed in RPA are 2-3 times 
larger than in the uncorrelated system (Reddy et al. 1999). This large 
suppression, which is similar to those encountered in the non-relativistic 
models, is due to repulsive forces in the spin-isospin channels. 

As mentioned earlier, RRPA is the self-consistent response of RMF 
models with vector and scalar interactions. It ensures thermodynamic 
consistency, i.e, it satisfies the compressibility sum rules discussed in 
§6.1.2. However, the cross sections are not greatly changed by interac- 
tions. In contrast, interactions in the spin-isospin channel do not affect 
the Hartree or mean-field ground state, but modify the neutrino opacity 
calculated in RPA. This underscores the need for detailed studies of the 
spin and spin-isospin response in nuclear models. To date, supernova 
simulations have employed equations of state (EOS) that are based on 
the mean-field approximation, including non-relativistic EOS's such as 
those due to Lattimer & Swesty (1991). Spin and spin-isospin suscepti- 
bility of these models remain largely unconstrained by either empirical 
data or input from microscopic calculations of nuclear matter. 

7. e+e - Annihilation 

Ignoring phase-space blocking of neutrinos in the final state and taking 
the relativistic limit (m e — > 0), the total electron-positron annihilation 
rate into neutrino-antineutrino pairs can be written in terms of the 
electron and positron phase-space densities (J-) (Dicus 1972): 

Qu ^ =Kl {^) (/L) / Jfe-Fe + (4-4 + +4-4 + )fce-de e+ , 

(118) 

where K t = (l/l8ir 4 )cer (C'J + C'l). Again, C' v = 1/2 + 2sin 2 #vK for 
electron types, C v = — 1/2 + 2 sin 2 9w for and v T types, and = 
(1/2) 2 . Rewriting eq. (118) in terms of the Fermi integral F n (rf), we 
obtain: 

Q Ue9e = Ki (kT) (-^) C^.) 6 [F 4 ( Ve )F 3 (- Ve ) + F 4 (- Ve )F 3 ( Ve )] , 



m e c 2 J \hc J 

(119) 

where r\ e = fi e /kT and 
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Integrating eq. (118), we obtain 



Q, 



< UePe ~ 9.7615 x 10 



,24 



kT 



f(r) e ) ergs cm 3 s 



i 



(121) 



MeV 



where 




(122) 



Q 



~ 4.1724 x 10 



,24 



kT 



f(r] e ) ergs cm 3 s 



4 



(123) 



MeV 



One can easily derive the spectrum of the total radiated neutrino 
energy (ex) by inserting a delta function (J 5(et — e e - — e e +)deT) into 
eq. (118). Recall that the total energy of the neutrinos in the final state 
is equal to the sum of the electron and positron energies in the initial 
state. Integrating first over e e + to annihilate the delta function and then 
over e e - to leave a function of et, one obtains: 



The numerical evalution of eq. (124) is straightforward. The average of 
ex is equal to: 



which near rj e ~ is ~ 8T and for r\ e >> 1 is ~ 4T(1 + r/ e /5). 

However, while the total energy loss rate (eq. 121) and the spectrum 
of et pose no great mathematical problems, the production spectrum of 
an individual neutrino is not so easily reduced to a simple integral or to 
an analytic expression. This is due primarily to the awkward integration 
of the angular phase space terms, while subject to the momentum con- 
servation delta function, and to the explicit dependence of the matrix 
elements on the electron/neutrino angles. From Dicus (1972), averaging 
over initial states and summing over final states, the matrix element for 
the e + e~ — > vv process in the m e = limit is: 



- A Y.\ M \ 2 = 16G2 ^ C 'v + C 'a) 2 P ■ 0* P' ■ Qi, + (C'v - C A ) 2 p • p' • qp] , 





(125) 



s 



(126) 
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where p and p' are the four-momenta of the electron and positron, re- 
spectively, and q v and qp are the four-momenta of the neutrino and 
antineutrino, respectively. Using the formalism of Bruenn (1985) and 
Fermi's Golden rule, expanding the production kernel in the traditional 
truncated Legendre series, performing the trivial angular integrals, tak- 
ing the non-trivial angular integrals from Bruenn (1985), and ignoring 
final-state neutrino blocking, we obtain for the single-neutrino source 
spectrum due to e + e~ annihilation: 

dQ 8vr 2 , f°° , 2 , , . 

where 

$> p (e u ,e p ) = — / de e -J r e -[e e -]J r e +[e 1/ + e D -e e -]H (e 1/ ,e i? e e -) , 
^ Jo 

(128) 

and 

Ho(e u ,Sp,e e ~) = (C v +C A ) 2 jQ(e u ,ep,e e -)+(C v —C A ) 2 jQ I (e u ,ep,e e -) . 

(129) 

The Jqs in eq. (129) come from the more obdurate angular integrals 
required by the dot products in eq. (126) and the momentum delta 
function and have the symmetry: 

Jq ( £ v> £ P> e e" ) = J^i^v-, £ e - ) • (130) 

From eqs. (127) and (129), we see that the differences between the 
spectra of the v e and neutrinos flow solely from their correspondingly 
different values of (C v + C A ) 2 and (C v — C A ) 2 . One can use 4-point 
Gauss-Legendre integration to calculate eq. (128) and 16-point Gauss- 
Laguerre integration to calculate eq. (127). 

At small rj e , the e + e~ annihilation spectra and total energy loss rates 
for the v e and v e neutrinos are similar, as are the average emitted v e 
and v e neutrino energies. However, as r\ e increases, both the total en- 
ergy radiated in v e neutrinos and the average D e energy start to lag the 
corresponding quantities for the v e neutrinos. This is true despite the 
fact that the total number of v e and v e neutrinos radiated is the same. 
If final-state blocking is ignored, (e.j)/T is a function of rj e alone, be- 
coming linear with i] e at high r/ e and one half of eq. (125) (~4.0) at 
low rj e . Note also that (s v )/T and (ep )/T are closer to one another 
than are {s Ve ) /T and (ep e )/T. The individual production spectra vary 
in peak strength, in peak energy, and in low-energy shape, but they are 
quite similar on the high-energy tail. Due to the parity-violating matrix 
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element for the e + e~ — ► vv process and the fact that r\ e is positive, the 
antineutrino spectra of all species are softer than the neutrino spectra. 
The pair sums of the integrals under these curves are given by eqs. (121) 
and (123). For rj e = 0, 50% of the pair energy emission of electron types 
is in v e neutrinos, but at rj e = 10 only 42% of this total energy is in u e 
neutrinos. However, at r\ e = 10, the v^ neutrinos still constitute 48.5% 
of the v^lvn pair emission. These differences reflect differences in the 
corresponding coupling constants C' v and C' A . 

8. UiUi Annihilation 

In the limit of high temperatures and ignoring electron phase space 
blocking, the Vi&i annihilation rate into e + e~ pairs can be written (Janka 
1991): 

Q v . p . = 4Kiir 4 (^2) J j ' JviJvMn + epj de n de Vi , 

(131) 

where J v is the zeroth moment of the radiation field, e u is the neutrino 
energy, Ki is defined as before {i.e., Ki = (l/187r 4 )c<7 (Cy + Cj 2 )), and 



3 

<J> '((^) ) (V*i),Pvi,Pi>i) = j 



(132) 

where the flux factor (n Ui ) = H v /J v and the Eddington factor p u = 
(p% ) = Pu/Jy Eq.(131) can be rewritten in terms of the invariant 
distribution functions T v : 

(133) 

Note that when the radiation field is isotropic ($' = 1) and when 
rj e = the total rate for e + e~ annihilation given in eq. (118) equals that 
for ViVi annihilation given in eq. (133), as expected. Buras et al. (2003a) 
have addressed the related and interesting process of v{Di — > VjVj. We 
refer to that paper for a discussion of the relevance and rates of this 
process. 

9. Nucleon— Nucleon Bremsstrahlung 

A production process for neutrino/anti-neutrino pairs that has re- 
cently received attention in the supernova context is neutral-current 
nucleon-nucleon bremsstrahlung (n\ + 712 — > 113 + n± + vv). It im- 
portance in the cooling of old neutron stars, for which the nucleons are 
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quite degenerate, has been recognized for years (Flowers 1975), but only 
in the last few years has it been studied for its potential importance 
in the quasi-degenerate to non-degenerate atmospheres of protoneutron 
stars and supernovae (Suzuki 1993; Hannestad & Raffelt 1998; Burrows 
et al. 2000; Thompson, Burrows, & Horvath 2000). Neutron-neutron, 
proton-proton, and neutron-proton bremsstrahlung are all important, 
with the latter the most important for symmetric matter. As a source of 
u e and v e neutrinos, nucleon-nucleon bremsstrahlung can not compete 
with the charged-current capture processes. However, for a range of tem- 
peratures and densities realized in supernova cores, it may compete with 
e + e~ annihilation as a source for v^, v T , and v T neutrinos ("f M "s). 
The major obstacles to obtaining accurate estimates of the emissivity of 
this process are our poor knowledge of the nucleon-nucleon potential, of 
the degree of suitability of the Born Approximation, and of the magni- 
tude of many-body effects (Hannestad & Rafflet 1998; Raffelt & Seckel 
1998; Brinkmann & Turner 1988). Since the nucleons in protoneutron 
star atmospheres are not degenerate, we present here a calculation of 
the total and differential emissivities of this process in that limit and 
assume a one-pion exchange (OPE) potential model to calculate the nu- 
clear matrix element. For the corresponding calculation for arbitrary 
nucleon degeneracy, the reader is referred to Thompson, Burrows, Sz 
Horvath (2000). The formalism we employ has been heavily influenced 
by those of Brinkman & Turner (1988) and Hannestad & Raffelt (1998), 
to which the reader is referred for details and further explanations. 

Our focus is on obtaining a useful single-neutrino final-state emission 
(source) spectrum, as well as a final-state pair energy spectrum and the 
total emission rate. For this, we start with Fermi's Golden Rule for the 
total rate per neutrino species: 

«- - < 2 ">7[nH] iSklSkry I>i*c>*~.. 

where 

E brems = ^ 2 (l-T 3 )(l-J^ 4 ), (134) 

<5 4 (P) is four-momentum conservation delta function, u> is the energy 
of the final-state neutrino pair, (uj u ,q u ) and (cop,q p ) are the energy and 
momentum of the neutrino and anti-neutrino, respectively, and pi is 
the momentum of nucleon i. Final-state neutrino and anti-neutrino 
blocking have been dropped. 

The necessary ingredients for the integration of eq. (134) are the 
matrix element for the interaction and a workable procedure for handling 
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the phase-space terms, constrained by the conservation laws. We follow 
Brinkman & Turner (1988) for both of these elements. In particular, we 
assume for the n + n^n + n + vv process that the matrix element is: 



where the 4 in the denominator accounts for the spin average for identical 
nucleons, Gf is the weak coupling constant, / (~ 1.0) is the pion- 
nucleon coupling constant, qa is the axial-vector coupling constant, the 
term in brackets is from the OPE propagator plus exchange and cross 
terms, k is the nucleon momentum transfer, and is the pion mass. 
In eq. (135), we have dropped q u • k terms from the weak part of the 
total matrix element. To further simplify the calculation, we set the 
"propagator" term equal to a constant (, a number of order unity, and 
absorb into £ all interaction ambiguities. 

Recently, Hanhart, Phillips, & Reddy (2001) have addressed these mo- 
mentum terms in the context of axion emission and v^a production in 
supernovae. In an effort to make contact with the approximation to the 
matrix element we present here, they plot Q as a function of average rel- 
ative thermal nucleon momentum (p; Phillips, private communication). 
The function peaks for ((p) between 150 - 200 MeV at C - 0.47. At 
p = 50 MeV C ^ 0.08 and at p = 500 MeV C ^ 0.27. We are most inter- 
ested in the region around the neutrinospheres, where the emergent 
spectrum might be most affected by nucleon-nucleon bremsstrahlung. 
Mass densities and temperatures in this region mig ht be 10 12 - 10 13 g 
cm~ 3 and 5 — 10 MeV, respectively. We estimate p in this regime to 
be ~ 175 MeV and take ( = 0.5 for all thermodynamical points. The 
constant A in eq. (135) remains. 

Inserting a j S(uj — uj u — ujp)du> by the neutrino phase space terms 
times ujuj v u}y/uj 2 and integrating over uj d yields: 



where again uj equals (lo u + ujp). If we integrate over oj u , we can derive 
the uj spectrum. A further integration over to will result in the total 
volumetric energy emission rate. If we delay such an integration, after 
the nucleon phase space sector has been reduced to a function of uj and 
if we multiply eq. (134) and/or eq. (136) by lo u /lo, an integration over to 




= A 



UJ u UJp 



(135) 




(136) 
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from <jj v to infinity will leave the emission spectrum for the single final- 
state neutrino. This is of central use in multi-energy group transport 
calculations and with this differential emissivity and Kirchhoff's Law 
(§2) we can derive an absorptive opacity. 

Whatever our final goal, we need to reduce the nucleon phase space 
integrals and to do this we use the coordinates and approach of Brinkman 
& Turner (1988). We define new momenta: p + = (pi + p2)/2, p_ = 
{pi — p2)/2, p3 C = ps — p + , and p^ c = p^ — p + , where nucleons 1 and 2 
are in the initial state. Useful direction cosines are 71 = p + -p-/\p + \\p-\ 
and 7 C = p + -pzcftp+WPScY Defining Uj = p 2 /2mT and using energy and 
momentum conservation, we can show that: 

d 3 pid 3 p2 = 8d 3 p + d 3 p- 

lo = 2T(u_ - u 3c ) 

u\,2 = u + + u- ± 2(u + u_) 1 / 2 7i 

u 3 ,4 = u + + u 3c ± 2(n + n 3c ) 1/2 7 c . (137) 

In the non-degenerate limit, the ^1^2(1 — ^3) (1 — ^4) term reduces 
to e 2y e~ 2( - u++u -\ where y is the nucleon degeneracy factor. Using eq. 
(137), we see that the quantity (u+ + u_) is independent of both 71 
and 7 C . This is a great simplification and makes the angle integrations 
trivial. Annihilating d 3 p^ with the momentum delta function in eq. 

(134), noting that pjdp = ( 2m p f v^J 2 dv,i, pairing the remaining energy 
delta function with u_, and integrating u + from to 00, we obtain: 

dQnb = gxirTL^ 7 ' 5 ^^ 

(138) 

The variable x over which we are integrating in eq. (138) is equal to 
2u3 C . That integral is analytic and yields: 

poo 

/ e~ x {x 2 + xuj/T) 1/2 dx = rje^K^r]) , (139) 
Jo 

where K\ is the standard modified Bessel function of imaginary argu- 
ment, related to the Hankel functions, and r/ = u/2T. Hence, the to 
spectrum is given by: 

^ oc e-" J l 2T uj b K 1 (uo/2T) . (140) 
dio 

It can easily be shown that [uS) = 4.364T. Integrating eq. (138) over 
uj and using the thermodynamic identity in the non-degenerate limit: 

/ 2ir \ 3 / 2 

e»=(^) ,„/2, (141) 



duo . 
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where n n is the density of neutrons (in this case) , we derive for the total 
neutron-neutron bremsstrahlung emissivity of a single neutrino pair: 

Q nb = 1.04 x l0 30 C(X nPl4 ) 2 (^) 5 - 5 ergs cm- 3 s" 1 , (142) 

where pu is the mass density in units of 10 14 gm cm~ 3 and X n is the 
neutron mass fraction. Interestingly, this is within 30% of the result in 
Suzuki (1993), even though he has substituted, without much justifica- 
tion, (1 + u/2T) for the integral in eq. (138) ([1 + (7T7//2) 1 / 2 ] is a better 
approximation). The proton-proton and neutron-proton processes can 
be handled similarly and the total bremsstrahlung rate is then obtained 
by substituting A 2 + X 2 + ?jX n X p for A 2 in eq. (142) (Brinkmann 
& Turner 1988). At X n = 0.7, X p = 0.3, p = 10 12 gm cm" 3 , and T 
= 10 MeV, and taking the ratio of augmented eq. (142) to eq. (123), 
we obtain the promising ratio of ~ 5(. Setting the correction factor £ 
equal to ~ 0.5 (Hanhart, Phillips, & Reddy 2001), we find that near and 
just deeper than the neutrinosphere, bremsstrahlung is larger than 
classical pair production. 

If in eq. (136) we do not integrate over uj u , but at the end of the 
calculation we integrate over uj from w„ to oo, after some manipulation 
we obtain the single neutrino emissivity spectrum: 



dQ' 



nb 



doj u 



2G (yf yi 1°° Y^Wfe - Vu) 2 dn (143) 



= 2C(§^ 3 f 6 -^(e - O l,2 di , (144) 

where r\ v = uj v /2T, C is the normalization constant equal to 3x5 ^ x11 
(= 0.564), and for the second expression we have used the integral rep- 
resentation of K\{rj) and reversed the order of integration. In eq. (144), 
Q n b is the emissivity for the pair. 

Eq.(144) is the approximate neutrino emission spectrum due to nucleon- 
nucleon bremsstrahlung. A useful fit to eq. (144), good to better than 
3% over the full range of important values of r/ v , is: 

dQ' nb „ 0.234Q n6 /c^\2.4 _ x 

dio u t It J ' { ' 

Thompson, Burrows, & Horvath (2000) should be consulted for a de- 
tailed discussion of nucleon-nucleon bremsstrahlung for arbitrary nu- 
cleon degeneracy. 
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10. Conclusion 

The processes that have been described in this paper are essential ele- 
ments of the neutrino-driven supernova explosion mechanism. Coupling 
these with radiation-hydrodynamics codes, an equation of state, beta- 
decay and electron capture microphysics, and nuclear rates, one explores 
the viability of various scenarios for the explosion of the cores of massive 
stars (Liebendorfer et al. 2001ab; Rampp & Janka 2000,2002). Recently, 
Thompson, Burrows, & Pinto (2003) have incorporated this neutrino mi- 
crophysics into simulations of ID (spherical) core collapse and have in- 
vestigated the effects on the dynamics, luminosities, and emergent spec- 
tra of weak magnetism/recoil, nucleon-nucleon bremsstrahlung, inelastic 
neutrino-electron scattering, and a host of the cross section corrections 
described above. The character of the spectra reflect the opacities and 
sources. In particular, the energy hardness hierarchy from v e (softer) to 
Vn (harder) neutrinos is clearly manifest, as is the distinction between 
the v e pre-breakout and post-breakout spectra. 

To date, none of the detailed ID simulations that have been per- 
formed explodes and it may be that multi-dimensional effects play a 
pivotal role in the explosion mechanism (Herant et al. 1994; Burrows, 
Hayes, & Fryxell 1995; Janka & Midler 1996; Fryer et al. 1999; Fryer & 
Warren 2002). Be that as it may, an understanding of neutrino-matter 
interactions remains central to unraveling one of the key mysteries of 
the nuclear universe in which we live. 
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Figure 1. The scattering kernel 7? out (e v , e' v , cos 8) for f e — electron scattering as a 
function of cos# for e v — 20 MeV and e' v — 2, 10, 16, and 24 MeV, at a representative 
thermodynamic point (T = 5 MeV, p = 10 13 g cm~ 3 , Y e = 0.4). 




Figure 2. For the same thermodynamic point as used for Fig. 1, the scattering 
kernel (R out , thick solid lines) for ^—electron scattering as a function of cos 9, for e v = 
20 MeV and e' v = 2, 10, and 16 MeV. Short dashed lines show the first-order Legendre 
series expansion approximation to R° ut , which is linear in cos#; R° ut ~ (l/2)<E>o + 
(3/2)<E>i cos 6. The long dashed line shows the improvement in going to second order 
in cos6» by taking R out ~ (l/2)$ + (3/2)$i cos 6 + (5/2)$ 2 (l/2)(3 cos 2 6 - 1). 
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Figure 3. The I = term in the Legendre expansion of the ^—electron scattering 
kernel, &o(sv, e'„) (eq. (37)), for T = 6 MeV and n e — 20 as a function of e' v for e„ = 5, 
10, 15, 20, 25, and 35 MeV. Note that for any e„, the neutrino is predominantly 
downscattered. The magnitude of <&q{e v ,e' v ) and sign of (w) are to be compared with 
those in Fig. 5 for ^-nucleon scattering. 




Figure 4- The scattering kernel R° ut (e v , e' v , cos 8) for v e — neutron scattering as a 
function of cos 9 for e„ = 20 MeV and e' v = 18, 19, 21, and 22 MeV, at a representative 
thermodynamic point (T = 5 MeV, p = 10 13 g cm~ 3 , V n = 0.5). Note that although 
the absolute value of the energy transfer (|e„ — e^|) is the same for both e' v = 19 
MeV and e' v = 21, the absolute value of R° ut (20, 19, cos 6) is greater than that of 
i? out (20, 21, cos#), reflecting the fact that at this temperature the incoming neutrino 
is more likely to downscatter than upscatter. 
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Figure 5. The I = term in the Legendre expansion of the neutrino-nucleon scat- 
tering kernel, $o(e*,e'„) (eq. (37)), for T = 6 MeV and function of 
e' v for e„ = 5, 10, 15, 20, 25, and 35 MeV. Note that for e v = 5 MeV the neutrino 
is predominantly upscattered, while for e„ = 35 MeV the neutrino is predominantly 
downscattered. The magnitude of $o(e„,£^) and sign of (w) are to be compared with 
those in Fig. 3. 




Figure 6. Dynamic structure function of a plasma of ions as a function of energy 
transfer lu (measured in units of the plasma frequency ui p ~ 0.3 MeV) and fixed 
momentum transfer \q\ = 2n/L ~ 6 MeV (left panel) and |g| = 6n/L ~ 18 MeV 
(right panel). Statistical errors for the results of the molecular dynamics simulations 
are also indicated. 
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(a) 
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(c) 



Figure 7. Feynman diagram representing the Schwinger-Dyson equations: (a) mean 
field (Hartree) propogator; (b) dressed vertex for coupling to an external current with 
in the RPA; and (c) RPA polarization tensor or current-current correlation function 
for the coupling to external currents. Dressed propagators are thick solid lines, bare 
propagators are thin lines, dashed lines are strong/electromagnetic interactions, wavy 
line is the external (weak) current, the filled circle is the bare weak vertex, and the 
filled square represents the dressed weak vertex. 




Figure 8. Logio of the Gamow-Teller structure function versus to/q for an incident 
neutrino energy of 20 MeV, energy transfers, u>, of 6, 8, and 10 MeV, two values of 
the effective mass (m* = [0.75m n , 1.0m n ]) and two values of the density (p = 3 x 10 14 
and 10 15 g cm -3 ). A temperature of 5 MeV and a Y e of 0.3 were used, as was the 
default vgt (= 4.5 x 10~ 5 ). (Figure taken from Burrows and Sawyer 1998.) 



